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1. Introduction 


at Background 


Because of the environment in which they operate, the 
seakeeping characteristics of a ship are of paramount impor- 
tance when assessing its overall performance. In the past 
this aspect of a ship's performance had to be judged by the 
results of model tests conducted at a point in time well into 
the preliminary design phase. While these tests provide 
results of good quality, they were not obtained until the 
pending design was quite firmly established. In fact, the 
results obtained by model tests had the characteristic of 
being just that, results, rather than an important input into 
the design cycle. The obvious desire then would be to have 
a tool capable of providing accurate predictions of seakeeping 
performance based on the data available in the conceptual 
design phase. These predictions could then be used to in- 
fluence the selection of hull form coefficients, etc. prior 
to the time when the hull form is actually being generated. 

In 1975 Professors T. Loukakis and C. Chryssostomidis 
published the "Seakeeping Standard Series for Cruiser-Stern 


{1} 


Ships". This paper corrolates the seakeeping behavior, 


as predicted by computer model, of the Extended Series 60 








hull forms and sets forth a method by which the performance 
of this type of ship may be predicted based on five para- 
meters: Froude number, F; ratio of significant wave height 
to ship length, S; beam/draft ratio, B/T; length/beam ratio, 
L/B; and block coefficient, CB. With this procedure a 
designer can predict the relative merits of various candidate 
designs at a very early stage. This represents a significant 
capability. 

As a result of the work represented in reference [1] 
there is considerable interest in generating a similar sea- 
keeping series for contemporary cruiser/destroyer type hull 
forms. In order to do this in the fashion of reference [1], 
a representative sample of the ship type must be analyzed by 
computer model and then the results corrolated. It was 
this need for sample hull forms that provided the motivation 


for this thesis. 


ft o2 Thesis Content 


There are two aspects of hull form generation 
addressed in this thesis: first, hull form modification and 
second, mathematical lines representation and fairing. The 
technique of hull form modification developed in chapter two 
is based on the work of H. Lackenby reported in reference [2]. 


The essence of this method is that the sectional area curve 


10 








of an existing ship is redrawn in a systematic fashion to 
produce a curve with the desired values of prismatic co- 
efficient, Cyr and longitudinal center of buoyancy, LCB. The 
sections are then shifted longitudinally to produce a modified . 
form with these characteristics. In applying this technique 
to destroyer type ships there were several anomalies encoun- 
tered which required that the method of Lackenby be further 
modified. These modifications, with the pertinent background 
are contained in chapter two. 

The other aspect of lines generation which is addressed 
in chapters three and four is mathematical lines representation 
and fairing. Although the fairness of a hull form is not 
critical to the seakeeping analysis it is an unavoidable 
subject when considering computer aided ship design. In 
these chapters the use of parametric cubic splines and least 
squares curve fitting are addressed. While the parametric 
splines are shown to provide the capability of representing 
virtually any type of line, the least squares fairing technique 
is limited to use with curves Pepresentarie by single valued 
functions. The algorithms are, however, capable of fairing 
lines with infinite slopes at the end points. 

It is anticipated that the tools developed in this 
thesis could be readily fused into a single computer program 


with the capability of modifying an existing ship form to 


shal 


obtain a faired design with the desired coefficients of form. 
When this is developed it will be possible to generate 
rapidly any number of designs for subsequent performance 
analysis. The implications of this are discussed in chapter 


six. 
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2. Method of Hull Form Modification 


Qe. Background 


During the design of all but the most trivial engin- 
eering systems, it is incumbent upon the engineer that he 
or she formulate a model of that system. Additionally, 
the designer must continually refine the model with each 
successive iterative cycle so that the results are of 
sufficient detail to be meaningful [3]. One such model 
used during ship design is a geometric description of the 
ship's hull form. The most traditional manner of providing 
this information is by way of the lines drawing. 

The ship's lines drawing, more frequently referred 
to as the ship's lines, is a set of three orthogonal views 
of the ship's hull depicting the lines of intersection of 
various planes with the hull form. When viewed in conjunc- 
tion with one another, they provide the capability to 
spatially locate any point on the moulded surface of the 
ship. Figure 2.1, taken from reference [4], is an example 
of a lines drawing for a "Mariner"-class, steel hull cargo 
vessel. 

While the lines drawing, prepared manually by the naval 
architect and draftsman, has been the older and more tradi- 


tional means of depicting the geometric properties of a ship, 
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the advent of the high speed computer has provided great 
impetus for defining the ship's form in a mathematical 
format [5]. It is interesting however, that a very 
successful attempt was made at representing ship's lines 
mathematically by Admiral David Taylor in the early 1900's 
[6]. This will be expanded upon a bit later. 

When first confronted with the job of creating the 
lines of a new ship, the naval architect seeks a means of 
quantifying the expected form of the vessel so that he may 
strive to create an "optimum" design. These optimizing 
criteria generally take the form of requirements and restric- 
tions placed on the various coefficients of form, i.e., 

Coe Ci LCB, LCF, etc. However, there might also be 
requirements placed on certain specific regions of the 

ship. An example of this could be the shape of the midships 
section for a cargo vessel or the stern configuration 
dictated by propeller and rudder selections. Nonetheless, 
when the naval architect completes his candidate design, 

the important product will be a faired set of ship's lines 
meeting all the optimizing criteria previously established. 

The above procedure is clearly long and involved. For 
this reason much effort has been expended to develop hull 
form modification techniques. The objective of these pro- 
cedures is to utilize an existing, successful hull design, 


Or parent form, as a basis and then to alter this form ina 


is 








systematic fashion. This modified form should have the 
desired characteristics, and, hopefully, require little 
additional refairing. It is this fairing procedure, des- 
cribed in chapter 4, which requires a large proportion of 
the designer's effort. The remainder of this chapter 
addresses the modification techniques themselves. 

One of the oldest and most widely used methods of hull 
form modification is illustrated in Figure 2.2. In essence, 
the sectional area curve of an existing ship is altered by 
some arbitrary or systematic method to produce a curve which 
satisfies some criteria of the designer, usually prismatic 
or block coefficient and longitudinal center of buoyancy. 
The offsets for the new design are then obtained by taking 
the section in the parent whose ordinate in the sectional 
area curve matches the ordinate of the derived curve. This 
is represented by the movement of section A at position x, 
in the parent to section A‘ at position Xa in the derived 
hull form. Hence, it is merely a longitudinal repositioning 
of the existing sections. This method works reasonably well 
as long as the shifts are of "moderate" amount and the 
designer is prepared to accept the resulting profile and 


waterlines without alteration. 
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Figure 2.2 





The above method lends itself quite well to design 
without the aid of computers or other automatic computational 
devices. However, in recent years there has been much work 
done in the area of hull form modification with the use of 
high speed digital computers. In virtually all cases where 
computers are used, an effort is made to represent the ship's 
contours or surface regions [5] in a mathematical format. 

It is for this reason that the "Taylor Standard Series" is of 
interest. It was Admiral David Taylor who, in the early 
1900's, generated one of the first successful hull series 
based on representing the sectional area curve and design 
waterlines by fifth degree polynomials [6]. 

Another procedure of hull form modification utilizes 
specific transformation functions to alter various regions 


or characteristics of the hull [7]. This form of 
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modification provides the user with much greater control over 
the specific ship form than the method of shifting sections 
longitudinally as previously described. An interesting 
description of this type of procedure may be found in 
reference [7]. 

It is, however, the method of longitudinally shifting 
sections which was selected for development in this thesis. 
The reasons for its selection are twofold. First, preliminary 
work with destroyer-type ships conducted at M.I.T. during the 
summer of 1977, indicated that the results of the modification 
were quite realistic and not plagued by gross unfairness. 
Secondly, the method was tractable and readily adopted to the 
peculiarities of destroyers. Those peculiarities being 
principly the fact that this type of ship has no parallel 
middle body and also that the section of maximum area, in 
most cases lies at a location other than midship. 

The specific method of modification used is that of 
Lackenby [2] as subsequently modified by Moor [8], and then 
again by this author. Briefly, the developments presented in 
reference [2] are highly general, permitting the designer to 
vary the value of prismatic coefficient, Cyr and the longi- 
tudinal center of buoyancy, LCB, of a very wide variety of 
ships. Included was the capability of altering, or retaining 
unchanged, the parallel middle bodies of ships so configured. 


However, one serious drawback was that the method left no 


18 





control over the design waterline, and while this line might 
turn out fair, the longitudinal center of floatation, LCF, 
merely ended up where it did. It was to this problem that 

- Moor [8] was concerned. By addressing himself to both the 
sectional area curve and the design waterline in the manner 
of Lackenby, and then coupling the two procedures, he was 
able to obtain a derived form having the desired values of 


Cc LCB and LCF. 


p’ Sw! 
At this point only one minor problem existed with the 
method as it stood. For ships with keelrise fore or aft it- 
was possible to obtain unwanted oscillations in the ship's 
centerline profile. To eliminate these oscillations, this 
author extended the logic of Moor to include the ship's 
profile. In so doing, the designer may be assured of a 
derived form having not only the four desired characteristics 
and coefficients previously mentioned, but also the desired 
profile. The only restricting requirement, other than the 
fact that the changes in Cc. and C. be "moderate", is that 
for the method to be mathematically rigorous the section 
of maximum beam, sectional area and local draft must be 
coincident. If this isn't the case a small (#1%) unpre- 


dictable error, based on the parent hull design and the 


desired changes is introduced. 


9 





All of these relationships are developed in full in 
| the following section. The only other alteration to the 
_ methods of Lackenby and Moor was that the procedure had to 
be capable of accommodating destroyer-type ships whose 
maximum sections fill other than at midship. This change 


is also included in the derivations that follow. 


20 





22 Development 


2.2.1 The "One Minus Prismatic" Variation 


As a means of introducing the method of longitudinal 
repositioning of sections, the traditional "one minus 
prismatic" is first developed. This procedure enables the 
designer to modify the fineness of a parent form by expanding 
(creating in ships without), or reducing the region of 
parallel middle body. It is convenient for this, and the 
following derivations to refer to Figure 2.3 and the 
following definitions. It should also be noted that the 
sectional area curve is normalized with respect to both the 


value of maximum area and length of the half body. 


Figure 2.3-SA. Curve 
(with parallel middiebody) 
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For the 


For the 


66 = 


ép = 





parent design: 

the prismatic coefficient of the half body. 
the fractional distance from midships of 

the centroids of the half body. 

the fractional parallel middle of the half body. 
the fractional distance of any transverse 
section from midships. 

normalized area of any transverse section at 
longitudinal position x. 

derived form: 

the required change in prismatic coefficient 
of the half body. 

the resulting change in parallel middle body. 
the necessary longitudinal shift of the 
section at x required to generate the required 
change in prismatic coefficient. 

the fractional distance from midships of the 
centroid of the added "sliver" of area 


represented by 6¢. 


In Figure 2.3 it should be recognized that AB'C is the 


curve of the derived form and curve ABC is that of the 


parent. In accordance with the method of the "one minus 


prismatic" the new location of the transverse sections is 


defined by the following equations. 


Ze 





1-(x+6x) _ 1-(o+66) 
1=$ 


1l-x 


hills: 
ye 


6x 
i. 


§ 
ts (1-x) (2.1) 


6x 


The area BCD is seen to be 1 - ¢ and the above modification 


l - (+56) 
1-6 r 
The new area B'CD is therefore 1 -— (¢+66), demonstrating 


simply reduces it by the factor 


that the method generates the desired prismatic coefficient 
of ¢ + So. 
There is however a concomitant change in a parallel 


middle body found by solving for x at x = p, i.e., 


ép = S#. (1-p) 


Sp_ = So (2.2) 


Therefore the resulting change in prismatic coefficient is 
obtained by altering the length of parallel middle body and 
then proportionally expanding or contracting the entrance 
and run. Because of this procedure the method has the 


following disadvantages: 
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l. There is no control over the length of the 
parallel middle body, i.e., ¢ and p, cannot 
be varied independently. 

2. The procedure cannot be applied to reduce the 
fullness of a ship having no parallel middle 
body. 

3. Conversely, a ship cannot be increased in 
fullness without introducing parallel middle 


body. 





4. The prismatic coefficient of the entrance or 
run cannot be altered... . 

5. The region where fullness is added cannot be 
controlled. That is, the maximum changes in 
fullness take place at the shoulders of the 
curve, i.e., point B. 

It is because of these numerous, severe restrictions that 
Lackenby sought to develop a more general technique of 


modification. 
2.2.2 Varying the Fullness of an Entrance or Run not 
Associated with Parallel Middle Body. 


In reference [2], Lackenby concerned himself with 
providing a means by which to modify both Cpe LCB and the 


length of parallel middle body in a controllable manner. 
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While many of these relationships are of importance, only 
those which apply more specifically to destroyer-type hull 


forms will be pursued in any detail. However, for the 


readers'convenience, the most generalized case of Lackenby's 


formulas is included as equations (2.10) through (2.13) in the 
last part of this section. For the in depth derivations the 


reader is referred to the original paper. Nevertheless, the 


following derivation is the foundation upon which all the 


, 


subsequent relationships are based. 


In referring to figure 2.4 the various quantities have 


a meaning identical to those of the previous section. The 


only additional term requiring definition is k, defined 
mathematically as follows: 
8 


39 oe x” dy 


we 


The only other difference between figures 2.3 and 2.4 is that 
figure 2.4 represents a hull form not having parallel middle 
body, i.e., p = 0, and as a consequence the length of the 
entrance and run equals that of the half length of the ship. 

In referring to figure 2.4, it is recognized that in 
order to preserve the form of the parent at both the end of ° 
the ship, (x = 1) and the middle of the ship, (x = 0) an 


equation for 6x of the following form would suffice: 
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ox ="ex (1 = x) 


where c is an as yet to be determined constant. It may be 


seen in Appendix A that the relationship for 6x is: 


Be = Pe id c= x) (2.3) 
$ (1-2x) 


Figure 2.4-S.A.Curve 





Clearly, the equation for 6x is of the form of a second degree 
polynomial (parabola) whose values are 0 at x = 0 and x = 1, 

as desired. This relationship also shows that the amount by 
which any section in the parent is shifted is a function solely 
of the unchanged longitudinal position x and some as yet unknown 


value 64. 
26 





area under the sectional 


axis (say x = 0). Since 





to the entrance and run, 


| 
'6¢'s of these respective 


: 
f 


At this point we must turn our attention and consider 
both the entrance and run concurrently if we are to lend some 
significance to the quantity 5¢. If we desire to specify both 


Cy and LCB, we are in essence placing a requirement on the 


area curve and its moment about some 
equation (2.3) applies independently 
it should be possible to select the 


regions such that when taken together, 


the ship has the values of C_ and LCB desired. 


P 


At this point we introduce the following quantities: 


Z = the distance of the parent ship's LCB 


from midships, normalized by the half 


length, 


(positive forward). 


6z = the required shift in LCB to obtain that 


required for the daughter hull form 


(positive forward). 


Prime (') 


Subscripts 


- denotes derived forms. 


- e = entrance or forward half-body. 


r run or after half-body. 


t a property describing the 


entire ship. 
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im 
Figure 2.5-SA.Curve 





Referring to figure 2.5 above, we may interpret the require- 


ments on C and LCB mathematically as follows: 


P eo = 2 
dL = 6, + 66, = (6, +O.) + (80, + 66.) (2.4) 


z= 2+ 62 
Summing moments of areas 
z’ ($, + 6o,) = Xabe a hoo, = [x,o, im h69,] 
28 = or (ko, + hobo, - [eo, + A, 66,1) (2.5) 


Oe 
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If equations (2.4) and (2.5) are rearranged, the expressions 


for 56. and So. may be obtained. 


66, = asa {66,(h,+2) + 682(o, + 66,)} (2.6) 
2 —_ - 
aoe = hah) {6o,(h,-Z) - 62(6, + 66,)} (2.7) 


At this point the only variables which were not previously 
defined are h, and h- The exact expression for these 


variables is, with the appropriate subscript: 


2%-3k7 , 6¢ (R= 3k* + 2x), 


h= 
1-2 (1 = 2x) 


(2.8) 
While equation (2.8) contains 6, the very thing for which it 
is being used to calculate, it has been stated [2] that the 
leading term along provides a very good approximation to h 
for "moderate" values of 6é¢, i.e., 

= 2 

h= eae Oke (2.9) 

1 - 2x 
Should it be desired to calculate 6¢ using equation (2.8), the 
solution will prove to be a quadratic which, while unwieldy, 
is certainly not unsolvable. The derivation of h may be found 


in Appendix A with the value of r defined as follows: 
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While the above equations with *he derivations in the 


| 
, 
! 


| Appendix illustrate the underlying theory, the following 


expressions represent the most general form of Lackenby's 


work. 


é.=5 : 5 {2[5o, (B_+z) +62 (9, +66.) ]+C,6D,-C, 5D, } 
e es 
(2.10) 
= a Zz Zz - 
6o.. = B+ 5, {2[6o, (B,-z)-5z (¢,-66,)1-C, dp, C,5P,} 


(2411) 
In the following expressions the items refer to the entrance 


or run as appropriate: 


al 6p , (x-p) pc. ¢ (1-9) 
éx = (1 - x) i + LEP) [69-6p tissyl} (2.12) 


The practical limits on 6b. and 5. are: 


ép(1-9) +5ati - $2] 
66 = . i-p (2.13) 


A, B and C are calculated as follows: 


A = 9(1-2x) - p(1-$) (2.14) 

HS $ {2x Bea =) olen (2.15) 
ete Bey Sos 

C = pes{B(1-9) - ¢(1-2%) } (2.16) 
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It should be recognized that in the above equations the 


necessary section shift, 6x, is a function of the new values 
of Cp and LCB, and the properties of the original parent hull 


form. 


2.2.3 Modification of Lackenby's Method to Accommodate Hull 


Forms with Maximum Sections not at Midships. 


It should be realized that in all of the preceding 
developments it was assumed that the section of maximum area 
fell at midships. While this is certainly the case for a 
large class of vessels, it is virtually never true for contem- 
porary cruisers and destroyers. It is for this reason that a 


new set of equations was sought while still adhering to the 


basic philosophy of longitudinally shifting sections. 


. Se See eee 


Le 
Figure 2.6-SA.Curve 
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Referring to figure 2.6, the definition of terms is, 
Once again, consistant with the preceding section. There are, 
however, two very important changes. First, where z and 6z 
“were originally normalized by the ships half length, they are 
now normalized by twice that length or the length between 
perpendiculars, L. Second, x, the local longitudinal position 
of any section is normalized by the appropriate length of 
entrance or run L, or L. respectively. Also all values of 
x and z take as their origin the station of maximum sectional 
area.- 


The basic relationship for 6x is still of the form 


ro a pe: = i 
6x = ex(1l-x) or 6x = ery) x(l-x), the same as equation 
(2.3) previously. However, equations (2.4) and (2.5) now 
become: 
weg 
o, ST {Le (o, + 6.) + L.($, + 5o)} (2.37) 
Pes yo Tee x +56 h ene x +66 h_)} 
4 e’'e”e ee Per oe rr 
t (2.18) 


These two equations are solved simultaneously for 5$ and 


i in Appendix B, the results of which are listed below: 


mm 1 Dime 2, ee ee 2 rn, 
te ‘ Lh +b_L bh. eyed Ute a™ edie Lh, (Leo thy oy 1g )} 
e e e rr (2.19) 
= J 2 Dy eee 
6 = = - ' - = ' } 
>) ae Te {Ll .X, Lv? x, z'L oS Toa eC? atlas 1? t) 
ool ree 


(2.20) 
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It was these two equations along with equation (2.3) which 
proved to give very satisfactory results for several sample 


calculations. 
2.2.4 A Method by which Constraints may be Placed on the 
Design Waterline 


In the preceding development the designer had no 


control over the shape of the design waterline. Because the 


‘longitudinal center of flotation LCF, was felt to be an 
important parameter in determining a ship's performance in 

a seaway, Moor [8] further developed the method of Lackenby 
to include control over the design waterline. It was with 
his revised method that Moor and his colleagues developed 
four distinctly different models with only their midships 
section identical. They were thus able to cut these four 
models in half to generate sixteen uniquely different hull 
forms. 

An interesting side light of this experiment was that 
the parent form used was that of a fast twin-screw currently 
in service whose maximum section was abaft midships. They, 
therefore, had to first swing the original area curve to 
place the maximum section at midships and then proceed with 
the new method of modification. Although having the maximum 


section at amidship may have proven to be more tractable for 
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the creation of the models, it was not necessary for the 
application of Lackenby's method. This fact was demonstrated 
in the previous section. 


The essence of Moor's method is that the sectional 


area curve and the design waterline are both altered in the 


sense of Lackenby to produce the desired values of Cor LCB, 
Cy and LCF. At this point a new factor is introduced; the 
ratio between the sectional area ordinate and the design 
waterline ordinate is calculated for both the parent and 
derived hull forms. These ratios are then plotted as a 
function of ship length and it is from this curve that the 
longitudinal shift of sections is determined. Figures 2.7 
and 2.8 illustrate the sectional area and design waterline 
curves and the area/waterline ratio curve respectively. 
Referring to figure 2.8, to obtain the offsets for a 
particular section R 


d 
parent is used as a basis. The reason for selecting station 


in the derived form, section a in the 


d 
the same value of area/waterline ratio. The offsets of 


e in the parent is that it is the closest section to R, with 


section RS are then multiplied by the ratio of the beam 
coefficient in the derived form at section R, to that of the 
parent at section Ro: 1 .6a> B,/8.,: These values may be seen 
in figure 2.7. Additionally, if the maximum beam of the 
derived form is different from that of the parent, the 


offsets of section R. also have to be multiplied by the ratio 
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of the maximum beam in the derived form to the maximum beam 
/d 


for any offset ba in the derived form is: 


in the parent, i.e., b Therefore, the equation 


dmax’ pmax“ 


B. b 
ba = b = ae (2.22) 
P “p> “pmax 


It can be seen in figure 2.8 that there are regions of 
ambiguity. Such is the case where the derived curve lies 
below a minimum in the parent curve. It has been this 
author's experience confirming that in reference [8], that 
the regions which cannot be explicitly be defined may be 
faired in after defining the sections on either side. 

The one remaining undesirable characteristic occurs 
in regions where there is some form of keel rise, i.e., the 
fore foot or skeg region. In these areas, if the draft of 
the parent is proportionally altered to equal that of the 
derived form, there is a concomitant and undesirable change 
in the area of the section. It is to this matter which the 


next section is addressed. 
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2.2.5 A Method by which Constraints may be Placed on the 


Ship's Profile 


Since Moor's method proved capable of constraining 
both the sectional area curve and the design waterline, it 
was decided to extend the method to include the centerline 
profile of the ship. The actual mechanics require only the 
introduction of a local draft coefficient, (local draft/ 
maximum draft) into the denominator of ihe avenypeam ratio. 
This new ratio, (Area/Beam/Draft), is graphed and the sec- 
tional shifts determined from this graph. In determining 
the new offsets not only are the offsets of the parent 
modified transversely as described in the previous section, 
they are also altered in the vertical sense. This altera- 
tion is accomplished by using the water plane as a reference 
and moving the waterlines below a distance proportional to 
the ratio of the derived form draft coefficient/parent draft 
coefficient. Also if there is a difference in the maximum 
draft of the derived form and parent, the waterlines are 
altered by this ratio as appropriate. 

As was mentioned in the background section of this 
chapter, section 2.1, for this modification technique to 
be mathematically rigorous the sections of maximum area, 
maximum design waterline breadth and maximum draft must be 


coincident. If this is not the case, the actual areas of 
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the sections generated will be consistantly different by a 
very small amount from what is desired. From the few 
examples this author has worked, it is estimated the 
difference in the value of Cp obtained and that desired is 
on the order of 1%. The explanation of this is seen in 


Appendix C. 
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3. Mathematical Representation of the Lines of a Ship 


S..1 Background 


It was established in section 2.1 that before the 
naval architect attempts to actually draw the lines of a new 
ship, he must have a "firm" description of the new design as 
represented by the various coefficients and curves of form. 
Examples of these, as cited previously, are: the sectional 
area curve, and hence c and LCB, the design waterline curve, 
(C, and LCF), the principle dimensions and perhaps specific 
information about the geometry of the midship section or 
stern region. These characteristics should represent what the 
designer feels is the “optimum” solution to his set of re- 
quirements. The naval architect now has to create one or 
more, of a possible infinity of, design candidates which ful- 
fill his descriptive coefficients. 

The traditional method for drawing the various lines 
of the ship is with the use of long, continuous strips of 
wood, metal, or more recently plastic, held in the desired 
position by weights. These tools are called splines and 
ducks respectively. The curves produced by this method 
were continuous but often times contained unwanted waviness. 


Removing these unwanted undulations, while preserving the 
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desired character of the line, is a process referred to as 
fairing. This topic, and the implications of placing a 

Be matical interpretation on it, are discussed in the 

next chapter. Not only did the naval architect have to 
generate smooth curves which pleased him visually, there also 
had to be a consistancy in location of the surface points 
when observed from the different views. This is sometimes 
referred to as cross-fairing and is also addressed in the 
final chapter. It is the fairing, and cross-fairing, which 
represents a very large part of the manual design effort. 

It was recognized long ago, that if the ship design 
process was to be automated to any degree, a technique to 
represent the lines of the ship mathematically would have 
to be developed. This is especially true today where much of 
the work is to be done by high speed digital computers. Not 
only must the designer/programmer provide the mathematical 
algorithms for representing the ship's lines, he must also 
provide the logic necessary for the computer to duplicate the 
heretofore trial and error methods of the draftsman. The 
alternative to programming the logic however, would be to 
give the system an interactive man-machine interface at the 
decision points. It is, however, the mathematical repre- 
sentation of these ships' lines to which this chapter is 


devoted. 
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B22 Development 


3.2.1 Derivation of the Spline Cubic Equation using a 


Variational Approach 


There are essentially two different methods by which 
one may arrive at a mathematical representation for ships 
lines. 

1. Select some mathematical function with several 
unspecified parameters whose values may be deter- 
mined by some accuracy criteria and boundary 
conditions. Typical of this approach is the use 
of a polynomial and a least squares fit criterion. 

2. Choose some smoothness and closeness of fit 
criteria such that, when taken together with the 
boundary condition, the function and parameters are 
determined. 

It is this second method, based on a variational smoothness 
criterion, that will be developed in this chapter. 

In general these variational methods involve the 
Minimization of the integral of some linear combination of 
the squares of the various derivatives of the function 
sought. In the case where the equation of the flexible 
spline is sought, the "smoothness" criterion is taken to be 


the minimization of the strain energy in the spline. 
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Mathematically this may be represented by [9]: 


b 


fL yyy! sy", eeey 
a 


,x) dx = min (3.5) 


where, for the spline equation (3.1) becomes: 


f ck*ds (3.2) 


Ss 


s = the total path length 


flexural rigidity of the beam 


Q 
i} 


nn 
ll 


curvature defined mathematically as: 


1-y 


ds = elemental arc length 


sey. + y" dx 


If these values are substituted into equation (3.2) the 


smoothness criteria becomes: 


x ne 
teats. — 575 dx = min (3.3) 
xy (1t+y'™) 


To complete the variational problem one must also 
consider a “closeness of fit" criterion which takes the 


following form: 
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Plysy'ey"eeeey em ,£,%) dx (3.4) 


a 
N 
94.0 


It may be seen in most any text on the calculus of variations, 
e.g.-, [10] that the criteria for smoothness and closeness 

of fit may be combined by the introduction of the unknown 
Lagrange multiplier A{[10]. The results of combining 
equations (3.1) and (3.4) into a single variational 


problem is: 


- 


p ' " (n) 
6 f {Lly,y 1Y peeelY 7x) + 
a 


AP (y,y',y",+.ey@™ ,£,x)} dx = 0 (3.5) 


A necessary condition for the integral in (3.5) to be 


stationary is: 


Q(LHAF) _d_ Q(LAaF)) , a? a (Ly F) 


ty 7 ax 
na@™ 9 (D+F) 
ett) an te } = 0 (3.6) 
dx ay 
with boundary conditions: 
2°. 2 
a a: ae) a ey ~ OL ou do. 0 
anes ey” eeteoye 515 ay" Gade. by™™ " G2, iv 
x dx” oy 
aL d 9L (n- 1), 
eee) Sy’ + ee ee er aye, o. y 0 
gytt FX ay2V ta) - ~ (3.7) 
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Equation (3.6) is known as the Euler differential equation 
for the variational problem presented in equation (3.5). The 
usual procedure for solving this system of equations is to 
first solve the Euler equation (3.6) in conjunction with the 
boundary conditions, expressed in equation (3.7). This 

solution results in an equation of the form y = £(A,x). This 
equation is then substituted into the "closeness of fit", or 
accuracy criterion of equation (3.4), to determine the value, 

or values of i. 

Another aspect of the mechanics of this procedure is 
revealed when the accuracy criterion requires that the 
resulting relationship for y = f(x) pass through a discrete 
set of data points. Such is the case for a "colocative spline", 
or a spline made to pass exactly through discrete data points. 
In this instance the integral in equation (3.4) would become 
a summation. However, in order to preserve the consistancy 
and similarity of working with integrals in both portions of 
equation (3.5), the Dirac delta function may be introduced 
anco F. 

Redirecting attention to equation (3.3), it will be 
noted that due to the complexity of the integrand, the 
differential element of strain energy, the result of equation 
(3.5) will not be closed form. In order to simplify the above 
2)5/2 


integrand it is assumed that the demonimator, (l+y' is 
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approximately 1. Or mie is very small. While this may not 
be the actual case, if the value of y' is linearized as being 
the value of the slope of the chord between two successive 
data points, the integral in (3.3) might be thought of as 
follows: 


2 


™ W(x) y"* dx = min (3.8) 


x 
I* = f 
aM 
If we were to take some mean value for W(x) for the entire 
domain of x, the minimization would be similar to the 
minimization of: 


- 2 
gu oy"* ax (3.9) 


xy | 
It is also this simplification which will allow us to calcul- 
ate a closed form solution to y = f(x). 

One other simplification which will be made, without 
harm to generality, is to assume that xe[0,1]. That is 
O <x) < Xo < eee SX) <x, <1. At this point we must 
actually define the accuracy criterion. Assuming that the 


curve passes exactly through the data points, we may say: 


y (x; ) = £5 i= 1 em (3.10) 
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Translating the conditions of (3.9) and (3.10) into the 


single variational problem of the form of (3.5) we have: 


1 2 m 
6 f fy"“+2 2 A.6(x-x.) (f-y)} = 0 (3.11) 
0 j=l ? y 


It should be recognized that the first delta (to the left of 
the integral sign) is the symbol for the variation, while the 
second is the Dirac delta function. f is any "candidate" 
function passing through the data points (x5, y(x5)) and 

the A's are to be determined from the accuracy criteria 

of equation (3.10). 


The resulting Euler equation (3.6) is: 


(eee Sez) 
“5 A, d(x x5) 0 ( 


The solution of this equation generates a function of the 
following form: 


a 


m 
y(A,x) = pet [> + ax? + Bx? 
Di 


+ Ax + Bx 


5 dy [xox 


i 


ees LD (3213) 


Therefore the value of y and the integration constants A, B, 


C, and D are all linear functions of the A'S. 
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The boundary equation (3.7) for this specific problem 
turns out to be of the following form: 


i 
[ey OV tay OY * ) = 0 (3.14) 
0 


It must also be noted that (3.13) is valid for 0 = oe 

While it is not intended to determine equation (3.13) 

for every possible situation, suffice it to say that the 
values of the A, 8 and the constants A, B, C, and D are 
uniquely determined by the coordinate points and the end 
conditions of the curve [9]. The essential points to be 
made are: 

1. For the criteria used, the equation for the 
spline as obtained by the variational approach 
is of the form of a multi-coefficient third 
degree polynomial. 

2. Where the curve is defined over the region (0, 1) 
by m data points, there is also a need for four 
additional pieces of information to satisfy, and 
fully solve equation (3.13). 

It will now be demonstrated that the form of equation 

(3.13) is also supported by the theory developed for the small 
deflection of elastic beams. As it turns out, the simplifying 


assumptions made in the preceding derivation are exactly that 
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which will be made in the small deflection theory. Neverthe- 


less, the consistancy of results lends much reassurance. 


3.2.2. Deflections due to Bending of a Simple Elastic Beam 


In virtually any undergraduate strength of materials 
course the subject of beam deflections may be presented in 
several different ways [11]. However, it will be the 


method of multiple integration which will be developed here. 


Figure 3.1-Strained Beam 
Element (Pure Bending) 





Referring to figure 3.1 above, it may be said that 
the element of the beam deforms about the neutral axis and 


that as a result, transverse plane sections remain plane after 
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deformation. This results in an elongation of those fibers 
outside of the neutral axis and a compression of those fibers 
inside. Also, the amount of distortion is proportional to 
the distance from the neutral axis. This being the case, 


the following equations hold: 


21a 
| 
tyo 
| 
m 
1 
ta} o> 
| 


therefore: 


(3.15) 


Te 


Le 
) 


In the above equations the following definitions apply: 


M bending moment 

E = the stiffness or Young's modulus 

I = the moment of inertia of the cross section 

p = the radius of curvature of the beam, measured to 
the neutral axis 

EI = the "flexural stiffness" or "rigidity" 


The mathematical expression for the radius of cruvature, or 


more traditionally the curvature, K, is defined as follows [12]: 
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Dir 


J es 
dy tea72 
[1+ (2X) 7} 


(3.16) 


It is to this equation which the simplification is applied. 
For actual beams it is assumed that the value of =2 is very 
small, hence, (GX)? <<l. If this is the case, the expression 


for curvature becomes: 
= =x=4 (3217) 


Therefore, when equations (3.15) and (3.17) are combined, 


the results are as follows: 


2 
E124 = M(x) (3.18) 


dx 


Here, M(x) is meant to indicate that the bending moment is a 
function of x. 

At this point we must address ourselves to the 
equation for the bending moment in a beam. Specifically, 
we will look at the results of loading a uniform elastic 
beam with concentrated point loads; remembering that this 
case most closely approximates the naval architect's ducks 


and splines. 


Sa! 





Figure 3.2 - Discretely 
Loaded Beam 





It is advantagious, at this time, to introduce the 
concept of the singularity function defined as follows: 
< 
o; x = Xs 


<x-x,>" = (3.19) 


n 
X~X. X > X, 
( i? BE 


With the aid of the singularity function, and in reference to 
figure 3.2, the bending moment equation obtained from the 


application of concentrated point loads is: 
M(x) = Pox -b Pi<x - x,> + P.<x - Xo> + ia. 


P <x - Xn? (3.20) 
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It should be obvious that when equation (3.20) is 
substituted into equation (3.18) the result may be twice 
integrated to obtain an equation for y = £(x) which is of 
the following form: 

3 3 3 


{P x + P,<x - X> + P<x - X5> 


air 


Ely (x) = A + Bx #+ 


3 


The above equation is of a form very much the same as equation 
(3.13). The essential difference is that in equation (3.13), 
the values of the end forces Ps and Ee were still unknowns 
requiring the statement of two conditions at each end of the 
beam. It should also be apparent that the resulting values 
of the A,'s are nothing more than 2EI times the forces 
required to keep the beam in equilibrium. Therefore, based 
on the results of this and the previous section, it will be 
accepted without further discussion that the third degree 
polynomial, or cubic, is an adequate model of the ducks and 
splines of the naval architect. Hence the term "spline 


cubic". 
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3.2.3 Piece-wise Continuous Cubic Polynomial Approximation 


Figure 3.3- Spline Curve 
Fitting Routine. 





Referring to figure 3.3 above, it is desired to 
approximate the curve y(x) by some series of cubic functions 


of the form: 
ga (x) = a. (x-x,) 2+b, (x-x,) 740, (x-x,) +d, (ane) 
) J e J J i 3 J 
where j represents the eis interval bounded by Xs and Xie1! 
and 1 < j Sn-1,n being the number of data points. In 


order for these segments to be continuous we impose the 


following conditions: 
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G5 (xj) = ¥ 


(1) j = Li2y660Nnk 


Sy-1 (Xp) = Yn 


(2) 95 (8544) * F541 (543) j=1,2,...n-2 
(4) G8(X449) = 84) (544) Jaye con=2 


From equation (3.22) it should be recognized that to fully 
describe the curve requires 4(n-1) unknown coefficients. 
However, the above equations provide 4n-2 conditions. We 
therefore require two more conditions. The obvious choice 
for these two additional constraints would be to specify the 
end conditions for the beam. Specifically, you would specify 
either g' or g" at the ends. 

For the sake of brevity the remainder of this deriva- 
tion will be abridged to include only the essential equations. 
For a complete and detailed description of this procedure the 
reader is referred to references [13] and [14]. 

Continuing with the derivation, the following 


definitions will prove useful. 


ha. = xX. = Xe (3.23) 
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D; = (Youa = y4)/h, (3.24) 


We may now relate the unknown coefficients in the following 


manner : 

ee is, = 9S2) (3.25) 

j oh. jtl1 j i 
2 
b, a (3.26) 
ny 

C5 = D5 Tog (23. * 5441) (3.27) 

Cra vee 3.28 
j Y5 ( ) 


Substituting these equations into condition (3) will generate 
a relationship between successive values of s. of the 
following form: 

Pee + 2th.) + ~ ane 

sh. (he h s 3) 


441) 8541 + 5542 


a eS Lp 2y<e oN=2 (3.29) 


Equation (3.29) will thus generate the following system of 


equations: 
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2(h, +h.) ae ie, 6(D,-D,)=s,h, 
hy 2(h +h) h, 8, 6(D3-D5) 
ee ee nee) nee 
ho-g 2¢h,_2th,_y)] $4 6(D,_4-D,_2)- 
| nate 
(3.30) 


It can be seen that, for any point x., 8, is the curvature 
a 
at that point. For the above system of equations if the 


curvature is known at the end points, i.e., s, and sn! the 


1 
curve will be completely defined. 

If instead of curvature the slope is specified at the 
end point, the above matrix will be modified erigheiye In 
this case, the value of the curvature at the end point will be 
unknown. For the situation where the beginning slope is 
specified the following changes will occur. From equation 


(3.27): 
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h 


1 a = 
D, eee (2s, + So) =o = ty 
or 2s,h) + Sohy = 6 (Dy - t,) (3.31) 


The effect on the matrix system will be to change the currently 
existing top row and then add another top row and left 


column as follows: 
2h, h, Reece Sy 6(D)-t,) 
hy 2(h, +ho) h, a ates So 6(D5-D,) 


If the end slope is specified a change of similar form takes 
place only adding a row to the bottom and a column to the 
right side as follows: 

|6(D 


2(h +h h s 


GP eee eee rt nomena. 


2h s 6(-D._ytt,) 


n-1 n 

The form of the above matrix is tridiagonal and lends 
itself to rapid solution by a recursive relationship [14]. 
This fact will save a significant amount of computational 


time when reduced to a computer algorithm. 
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3.2.4 The Rotating Spline [15] 


One particular disadvantage of the "piece-wise 
cubic" method developed in section 3.2.3 is that it will not 
provide a solution for curves having infinite slopes. It is 
for this reason that method of the "rotating spline" was 
developed. Referring to figure 3.4 it may be related that 
this procedure is merely a modification of the "piece-wise 


cubic" technique. 


Figure 3.4 - Rotating 
Spline Routine 





For this method of curve approximation a cubic 
polynomial is generated for each interval as before. However, 


in this case the coordinate system is redefined for each 
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interval, and the cubic equation is generated with respect 
to this local coordinate system. 

Appendix E contains the steps necessary to generate 
the computer algorithm. There is, however, one definate 
disadvantage to using a rotated coordinate system. In order 
to obtain an interpolated ordinate on the curve the 
following two parametric equations must be used. 

a= x, +(x “x, )t-(y;,)7y,) t(1-t) (a, (1-t) -b, ¢] 

(3.32) 


itl 


and 

y = Yst (ys 4) 7¥,) t+ (x5 17-x,) €(1-t) [a, (L-t)-b; t] 

(3.33) 

Both of these equations are third degree polynomials in the 
parameter t. To solve for some value y of the point (x,y) 
in the unrotated coordinate system, x is used in equation 
(3.32) to solve for t such that 0 < +51 and t is also real. 
This value of t is then used in equation (3.33) to calculate 
y- It should be pointed out that the quantities as and b; 
were determined previously as described in Appendix E. 

In summary, we have shown by two methods, variational 
calculus and simple beam theory, that the third degree or 
cubic polynomial provides a good representation of the thin 


elastic spline used in drawing ships lines. There was, 
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however, the disadvantage that the equations were not capable 
of representing curves having infinite slopes. For this 
reason the method of the rotating spline was introduced. 

This parametric method permits the representation of 
virtually any continuous curve, including those which are 


non singular. 
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4. Mathematical Fairing of Lines 


4.1 Background 


It has been stated previously that the lines fairing 
process can be the most time consuming aspect of the ship 
design cycle. For this very reason fairing becomes a prime 
candidate for automation. The difficulty, however, lies in 
the fact that obtaining a universally accepted mathematical 
definition of a faired line, or the fairing process itself, 
is a virtual impossibility. Perhaps the most general defini- 
tion, and one which would prove the least restrictive, is 
the following: 

A faired line is one which retains 
the desired "character" but eliminates 


any undesired waviness or fluctuations. 


It will be shown in the following sections that this 
result may be achieved by fitting, in a least squares sense, 
a third degree polynomial toa set of four or five data 
points. The number being dependent upon the desired boundary 
conditions. In addition to the similarity of the third 
degree polynomial to the form of an elastic spline, the 
polynomial also provides the capability of introducing a 


desired inflection point into a series of data points. It 
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also prohibits the introduction of multiple inflection points 
and undesired waviness, also an asset. Because of these 
characteristics and the excellent results demonstrated in 
reference [16], this "least squares" criteria was employed 


as the foundation of the fairing process. 
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4.2 Development 


4.2.1 The Least-Squares Criteria for Defining the Cubic Curve 


It was established in chapter three that the cubic 
polynomial would provide a "good" approximation of the shape 
of a spline used to construct the lines of a ship. What 
remains to be shown is how these polynomials are applied in 
order to generate the faired position of a set of data 


points. 


Curve II 


Curve I 


Figure 4.1- Least Squares Fit 





2 
x + a5x + a3Xx (4.1) 
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Referring to equation (4.1) above it should be 
recognized that, in order to uniquely specify the cubic 
polynomial in its general form, four independent pieces of 
information are required. The result of applying this 
information is to determine the values of Age ay ay and a3: 
This will produce a curve which exactly conforms to the given 
requirements. This is illustrated as curve I in figure 4.1 
where the curve is required to pass exactly through the first 
four data points. The disadvantage of using this type of 
curve is that, since it is required to pass exactly through 
the given data points, it is unable to modify their position. 
It is this alteration however that is necessary if the curve 
is to be "faired". 

As stated above, four pieces of information are 
required to uniquely specify the cubic polynomial. If, how- 
ever, we were to over specify the requirements of the curve 
and then demand that the solution satisfy these requirements 
in some "best possible" but not exact manner, we begin to 
get a feel for how fairing can be produced. Mathematically 
this can be stated as follows. 

Referring to curve II in figure 4.1 we will require 
that our resulting curve pass as close as possible to the 


five given data points. Usually this translates into a 
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mathematical form by requiring that the sum of the squares 
of the distances between the curve and the given data 


points should be minimized. 


ws 2 ae 
s=f ‘ [y; - (ay + aX; - AoXs + a,x) ] (4.2) 
L= 
or 
os. 9S) yoo a, 
day day da, 


These derivatives generate the following system of normal 


equations which can be solved for Ag cee 43> 
So 8) So 83 ap to 
aT So 3 S4 = Ss) 
a (4.3) 
a =3 24 25 a, =) 
83 S4 =e S¢ as ae 
Where Ss, and ty, are defined as follows: 
5 
Ss, = = me (4.4) 
i=l 
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y.x* (4.5) 


It can be seen that curve II in figure 4.1, while not 
passing exactly through the data points, passes "fairly 
close" and also displays a smooth and continuous character. 
It is this closeness of proximity, or minimization of the 
least squares difference, procedure that is the essence of 
the fairing criteria used in this thesis. 

The following sections will develop the equations for 
line segments whose end position and slope or just end 
position are fixed. First, however, a brief explaination 
of how these segments are applied to fair a complete line 


or set of data points. 


4.2.2 The Moving Strip Method 


In the procedure described above it was seen that a 
least squares spline was passed through five data points and 
the points on that line were then considered to be fair. In 
order to fair a complete set of initially unfair points 
toge 


k k+4° 
After passing a least squares spline through these five 


consider only five points at a time, i.e., P 


points obtain the faired position of point Pia? Then move 
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the strip one unit (k = k+l) and consider the next five 
points, PL to Pag! fairing P42? For each step we could use 
previously faired values for Py and Peed expecting our final 
solution to be obtained more rapidly. By walking this strip 
of five points through the entire set of data the faired 
position of each point may be obtained. This procedure may 
also be found in references [16, 17]. 

The following three sections develop the equations 
needed when considering data points whose boundary conditions 


are of the following type: 


1. Free end--the position and slope of the end is 
unspecified. 

2. Pinned end-~the end position is fixed but free 
to rotate. 


3. Clamped end--end position and slope are fixed. 


22.201 STRIP1L: Fairing an interval with free ends. 


This procedure is the same as that developed in 
section 4.2.1. It is this routine which is used to fair 
the center data point (Py 49) of five interior data points, 
i.e., PL. # Py and Pied # Po! where P and Pa are the first 
and last points respectively in the set of given data. This 
routine is also used to fair Pir Por Pa-l and 2 for the 
case where the ends are free to both rotate and translate. 
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and 


also 


Repeating equations 


So oF a0} 
24 55 =3 
29 =3 S4 
39 S4 5 
5 

k 
s, =f Soe 
k =) i 


. 2 
P(x) = ay + a,x + a5Xx +a 


(4.3) to (4.5) for convenience. 


= 29 al 
S4 i ei 
Ss a2 / Eo 
ae 233 — 


3 
3% 


(4.3) 


(4.4) 


(4.5) 


(4.1) 


Therefore the faired position of the second and third points 


in the five point strip are: 


Y2 


= P(x) and Y3 = P(x,) 


For the case of a free end point the first point becomes: 
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¥, = P(x) 


It would be fair to expect that the equation P(x) 
would be most representative of the actual curve at its 
interior regions where uncertainty about end conditions 
would have less effect. For this reason only the center 
point of a strip is recalculated as being faired, i.e., 


Prag a8 opposed to recalculating both PL and Peel’ 


4.2.2.2 STRIP2: Fairing an interval with pinned end. 


Figure 4.2 - Pinned End 
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In fitting the curve of equation (4.1) to the five 


data points in figure 4.2 above, we require that P(x) ) =y 


exactly. 


such that 


P(x) = Yi + a,x + anx? + a,x 


1 


If in our calculations we adjust the abscissas 


x 


0 we may simplify this equation to: 


3 


(4.6) 


Applying our least squares criteria to this we obtain the 


following normal equations: 


where 


Sk 


and 


4 
Ss 
Se 
xik 
1 
k 
= 1 
Cy TY) xy 
Sa 


Ti ert 
ee re, 
a, Be 
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(4.7) 


(4.8) 


(4.9) 





With the values of Ayr ay and ay computed from equation (4.7) 


we can calculate the faired position of points Po and P3: 


where 
: = _ 's = 
x 2 2 ae? = — a 


After these two points are determined STRIP1 can be applied 


to continue the fairing process 


4.2.2.3 STRIP3: Fairing an interval with clamped 


ends. 


Figure 4.3 — Clamped End 
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In order to fair an interval with both position and 
slope of the first point fixed we will consider only the 
first four points as shown in figure 4.3 above. In order 
to simplify the derivation we once again adjust the abscissas 


such that x) = 0. For this case equation (4.1) becomes: 


P(x) = y,; + ax + ajx* + a,x” (4.10) 


After applying the least squares fit criteria to the four 


data points the normal equations obtained are: 


S4 Ss 2 2 
a (4.11) 
ES 35 Ss a 
where 
4 
s, = xi (4.12) 
a ih 
i=2 
4 k 
= - ' — 
ty zo (y, qx; Yo) x, (4.13) 
1=2 
and 
's= = 
ant a Tl 
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1 


For the simple 2 X 2 system above a, and a, may be written 


as follows: 


t.s t.s 
a, = +t 33 (4.14) 
=) ars = ae 
t.s, - t.s 
a, a (4.15) 
2426 5 95 
where 
S4 Ss 
A = #0 
2s 56 


Using the values of as and a3 calculated the faired position 


of the second point may be readily determined as: 


Y2 


P(x) 


where 
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After fairing the second point STRIP1 can be applied to 
continue the fairing process, fairing point three and four 


the first time it is applied. 


4.2.2.4 TRANS1: Fairing the last points in a given 


sequence of data. 


As can be seen in the previous sections, the fairing 
procedures operate on a series of points with monotonically 
increasing abscissa and end conditions specified at Xi 


where KX) < Xo < Xz eee + At the time when the five point 


2 
fairing interval reaches the other end of the curve, i.e., 


PL = Pi-4! the following transformation must take place: 

t —_ — = t+ 

al xy Xn 0 Y Yn 
= = 's 

ta Xt Xn-1 Y2 Yn-1 

Xe] 4.0 xX. ey ee 
3 n n-2 : n-2 (4.16) 
t = = 

“ae nn a3 Y¥4 ~ Yp-3 
pS - i 

a5 xn *n-4 Y5 Yn-4 


This allows the fairing of the last three points as if they 
were the first three points. Once these three points are 
faired the reverse transformation is employed to place Yy 


into yn etc. 
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4.2.3 Fairing of Curves with Infinite Slopes 

In ship design it is not infrequent that lines are 
encountered which possess infinite slopes at one or both 
end points. Such is the case of a section through the bow 
of a ship equiped with a bulbous bow. Here if the offsets 
y are expressed as a function of z, an infinite slope will 
occur at the bottom of the bulb. While the parametric 
rotating spline of section 3.2.4 will accommodate such a 
form, the simple cubic polynomial of equation (4.1) will 
prove indeterminate for an interval containing an end point 


with infinite slope. 


Figure 4,4- Axis Rotation 
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Referring to figure 4.4 above, it can be seen that if 
the axis are rotated by some small anle 6, and the data 
points are redefined in this new coordinate system the 
fairing process may be carried out as normal. fThe following 


equations relate the coordinates in the two coordinate 


systems. 
x' = x cos® + y siné 
(4.17) 
y' = -x sin8 + y cos6 
x = x' cosé - y' sin6é 
(4.18) 
y = x' sind + y' cosé 


It is also assumed that for small values of 6, e.g., 10°: 


Ay = Ay', Ax = Ax' = 0 (4.19) 


In order to continue the fairing process the faired 


position of the first three points and the slope at the 


third point could be calculated in the rotated coordinate 
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system and then transformed back into the unrotated coordinate 
system. The faired position and slope of the third point 
could be used to continue the fairing in the unrotated plane. 


The following slope transformation is also helpful. 


dy _ ay' 
54 = tan {arc tan(3Z1) « 8} (4.20) 
Awe soe lL Other transformations. 


There are any number of transformations one could 
use to accommodate the problem of the infinite slope. One 


tried by this author was that of letting [18]: 


i 
x (1 - cos 6) 


“ 
It 


fA 
1A 


or 0 1 (4.21) 


@Q 
tt 


arcos (l - 2x) 


The curve is then plotted as a cubic in 6. This has the 
advantage of eliminating an infinite slope at x = 0; in fact 
dy/de = /r/2, where r is the radius of curvature of y = f(x) 
at x = 0. The disadvantage, as seen by this author, is that 
fairing will take place in the distorted y,® plane. Addition- 
ally, even though the resulting curves appear to be aesthet- 


ically pleasing, some apprehension exists regarding the use 
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of lines faired in the two coordinate planes. For this 
reason the author opted for fairing in a rotated coordinate 


system as opposed to one which was distorted. 
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5. Computer Algorithms 


ee Overview 


The system of subroutines developed in this thesis 
were designed to provide two distinct capabilities: (1) to 
provide a means of fairing a series of data points not pre- 
viously considered fair, and (2) to provide the capability of 
representing a series of data points by an analytical 
mathematical expression. This second feature would also 
provide a means by which slopes, curvatures, etc. could be 
determined by interpolation. The theory of these two pro- 
cedures was developed in chapters four and three respectively. 

The ultimate objective of these subroutines would be 
their utilization in a program to fair and draw an entire 
ship form. Because of this and the virtually infinite nature 
of the lines existing in a hull form, the program has to be 
capable of handling many line types. As an example of this, 
see figure 5.1, the programs require the following information 
as input. 

1. Independent variable coordinates. 
2. Dependent variable coordinates. 
3. Data point type. 


4. Number of data points. 
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5. Type of end conditions. 

6. End slopes if required. 

7. An indication as to whether the input data 

is fair as submitted. 
The only other pieces of information required for 
fairing are: 

1. TOL: Is a tolerance representing a limiting 
distance which any data point may be 
moved in the fairing process. 

.2. ACC: This number represents an accuracy 
which, if during the fairing process 
a point is not moved by more than this 
amount, it is considered to be ina 
faired position. 

3. LIMIT: This number sets a limit on the 


number of iterative cycles permitted 
in the fairing process. 


5.1.1 Specification of Point type 


The designation of point type is designed to be as 


consistant as possible with reference [19]. 


EOINT TYPE DEFINITION 


0 : Normal point at the beginning or 
in the interior of a continuous 
curved line segment. 


1 : Break point at the end of a continuous 
curved segment. At present this point 
is treated as if it were, pinned. The 
slope, while unspecified, is discon- 
tinuous. 
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Figure 5.1- Point Type Examples 
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5.1.2 Specification 


This point can be at the beginning, 
middle or end of a straight line 
segment. The slope is continuous at 
this point. This point must be 
specified where a curved segment 
joins a straight line segment since 
the point is considered to be a 
clamped end condition for the curved 
segment. 


Break point at the end of a straight 
line. The slope is discontinuous at 
this point. 


of End Conditions 


The end condition designation is made with a two digit 


real number of the following format, "B.E". Here B corres- 


ponds to the end conditions at the beginning of the line and 


E the end condition at the end of the line. 


END TYPE 


1 


ee 


DEFINITION 


Free end. The end point is free to 
both rotate and translate. 


Pinned end. The end point is free 
to rotate only, the position is fixed. 


Clamped end. The end is totally con- 
strained. It is free to neither 
rotate or translate. The slope must 
also be defined. 


Clamped end with infinite slope. The 
slope, whether + » is determined by 
the second data point, i.e., + © 


af Yo > por - = aft Yo < Yi° 
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5.1.3 Storage of Pertinent Line Data 


The information necessary to fully describe any line 
is stored in a 32 x 7 two-dimensional array. This array is 
labeled CRV in the subroutines and its elements have the 
following significance. At present the first thirty rows 
are for data point or interval information and the last two 
rows are for overall curve characteristics. This could be 
easily expanded to allow more input data. 

l. Colume 1, CRV(1,1) to CRV(30,1): 

Abscissa of the input data points. 


2. Colume 2, CRV(1,2) to CRV(30,2): 
Ordinates of the original data points. 


3. Column 3, CRV(1,3) to CRV(30,3): 
The faired ordinates of the data points. 


4. Column 4, CRV(1,4) to CRV(30,4): 
The point type, see section 5.1.1. 


5. Column 5, CRV(1,5) to CRV(29,5): 
The values of a; as defined in Appendix D. 


6. Column 6, CRV(1,6) to CRV(29,6): 
The values of be as defined in Appendix D. 


7. (Column 7, CRV(1,7) to CRV(30,7): 
The slope of the curve at the data points as 
defined in Appendix D. 
It should be noted that the elements of columns 5, 6 and 7 
are obtained as a result of the splinning option. The data 
in column 3 are obtained as a result of exercising the 


fairing option. 
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ELEMENT 3 DEFINITION 


CRV (31, 2) 


The number of data points. At present, 
6SCRV (31, 2) $30. 


CRV(31,3) : The slope at the beginning of the 
curve. Left blank if not specified. 


CRV(32,1) 3: An indication of the fairness of the 
curve. 
l1.=the data submitted is fair. 
2.=the data submitted is not fair. 

CRV(32,2) 3: End condition specification, see 
section 5.1.2. 

CRV(32,3) 3: The slope at the end of the curve. 


Left blank if not specified. 


The other elements of the CRV matrix are reserved for future 
use, e.g., in a full ship fairing program. 

There are three other aspects of the program which are 
of the utmost importance. As was seen in the development of 
chapter three, the mathematical curve representation, or 
splinning procedure is fully capable of accommodating multi- 
valued curves. The fairing option, however, requires that a 
curve be single valued over the domain of the independent 
variable. Therefore, while it is possible to fit a cubic 
curve to virtually any series of data points, care must be 
observed when exercising the fairing option. The second 
inviolable characteristic of the program is that the data 
points must be submitted in a monotonically sequential fashion, 


i.e., the points must not be submitted in a random fashion, 
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but rather as they are encountered while following the path 

of the curve. The last consideration is that in order to fair 
any curved line segment there must be at least six data 

points in the continuous curved region. This is true regard- 
less of the end conditions of the line as a whole or the end 


conditions for a line segment. 
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BS. Description of Subroutines 


A flow chart and subroutine listing may be found in 


Appendix F. 


5.2.1 Lines Fairing 


The subroutines included in this section are utilized 
to calculate the faired position of the given data points, 


i.e., column 3 of the CRV matrix. 


5.2.1.1 Subroutine PREFAR 


This subroutine takes the data in the CRV matrix 
and loads all the points on a continuous curve segment into 
three linear arrays; X( ), Y¥( ) and YORIG ( ) representing the 
abscissa, faired ordinate (the original ordinate for the first 
iteration) and the original ordinate respectively. This 
process is governed by the value of point type, CRV(I,4). 
With these arrays established PREFAR calls either FARCRV or 
FARLIN, depending on whether the curve segment begins with 
an infinite slope. 

Upon final return to this subroutine the values of the 
faired position of the data points will have been calculated 


and placed in column three of the CRV matrix. 
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5.2.1.2 Subroutine FARCRV 


This subroutine takes the data in the X, Y, YORIG 
array, for those line segments which have infinite slopes, 
rotates the coordinate axis 10°(1/18 radians) and then places 
the transformed points, equation (4.17), in an XPRIM, YPRIM 
and YOPRIM array. The subroutine then calls subroutine 
FARLIN to fair the first six data points in the rotated 
system. At this time subroutine SPLINE is called to determine 
the slope at the third point, also in the rotated system. 

The subroutine then completes fairing the remaining data 
points by matching the position and slope at the third point, 
in the unrotated system. That is, assuming point three to be 


clamped and beginning with STRIP3. 


5.2.1.3 Subroutine FARLIN 


This subroutine takes the points of a continuous 
curve segment and calls the various STRIP_ subroutines which 
actually compute the faired position of the points. FARLIN 
also calls subroutines FSTPTS and TRANS] to fair the first 


and last points in the sequence. 
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5.2.1.4 Subroutine FSTPTS 


This subroutine fairs the first three data points in 
a sequence of data points based on the end condition. STRIPI, 


2 or 3 are called as appropriate. 


§.2.1.5 Subroutine TRANSL 


This subroutine fairs the last three data points 
based on the end condition specified. Specifically it 


transforms the abscissa in accordance with equation (4.16). 


5.2.1.6 Subroutine STRIP1L, 2 or 3 


These subroutines are described in detail in sections 
4.2.2.1 to 4.2.2.3. They use, as arguments, the variables 
in the Xl, Yl and YO arrays. Additionally they require 
values of TOL and ACC which place limits on the amount which 
a point may be moved and the amount of movement 
which is considered to be negligable. For the case where the 
point would move by more than TOL from its original (unfaired) 


position, its movement is limited by the value of TOL. 


5.2.2 Lines Representation 


The methodology of representing a line by a parametric 
cubic equation was developed in chapter three. The actual 
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sequence in which the process is executed is described in 


the following sections. 


5.2.2.1 Subroutine PRESPL 


This subroutine examines the input data in the CRV 
matrix and places elements of continuous curved line segments 
into the X, Y and YORIG arrays. This assignment is based on 
data point type found in column four of the CRV matrix. 
Referring to figure 5.1A, the program would load the first 
eight points into X, Y and YORIG. Point nine, point type 5, 
would be used in conjunction with point eight to determine 
the slope of the curved segment ending at point eight. The 
program then calls SPLINE to carry out the actual curve 


fitting algorithm. 


5.2.2.2 Subroutine SPLINE 


This subroutine uses the data in X, Y and YORIG 
obtained from PRESPL and carries out the curve fitting 
algorithm presented in Appendix D. Although many intermediate 
terms are calculated, the only terms which are retained are 
ae b; and d.» these quantities are subsequently used to 


calculate interpolated values of the independent variable 


and slope at a point specified by the user. 
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5.2.2.3 Subroutine INTERP 


This subroutine determines the interval in which a 
desired value of a dependent variable is located. It then 
passes the coordinates of the surrounding points and the 
values of as and b for the interval to subroutine CALCY 
which calculates the value of the dependent variable and 


slope at the desired point. 


5.2.2.4 Subroutine CALCY 


This subroutine calls CALCT to obtain the value of the 
parametric variable T. With the value of T the interpolated 
value of the independent variable is determined. Since the 
dependent and independent variables are represented para- 
metrically, the slope of the curve is calculated by the 


chain rule as follows: 


sie 
" 
ale 
— 
A ou 
ry 


(5.1) 


Since the value of T is determined as being the root of a 
third degree polynomial, CALCY is designed to calculate the 
interpolated value of the independent variable and slope 
for up to three unique and real values of T. However, the 
subroutine is designed to print a warning that additional 
points are needed to specify the curve if T has more than 


one real value. 
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§.2.2.5 Subroutine CALCT 


This subroutine calculates the roots of the third 
degree parametric polynomial in T using the algorithm in 
reference [20]. This procedure is also presented in Appendix 
E for the readers' convenience. It should be realized, 
however, that only the real roots are calculated in the 
subroutine, the imaginary roots lack physical significance 
for the purpose of lines plotting. 

Appendix G contains an example of a data set that was 
first faired then splinned and then interpolated at points 
equal ta:ene~twentieth of the domain of the independent 
variable. Once again it should be emphasized that, in order 
to fair a curved segment, at least six data points must be 
defined in that segment, including the end points of the 


segment. 
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6. Conclusions and Recommendations 


6.1 Hull Form Modification 


When work on this thesis began, the initial goal 
was to develop a series of destroyer-like hull forms for 
future use in seakeeping analysis. Preliminary efforts, using 
the method of longitudinally shifting sections [2,8], while 
showing promise, indicated that additional work would be re- 
quired if the procedure was to apply accurately to destroyer 
type ships. Specifically, the method had to be adapted to 
ships whose maximum beam and section of maximum area did not 
lie at midships. These necessary changes were made success- 
fully and the method was also extended to provide control 
over the ship's celteftine profile. The primary motivation 
for this extension was to gain control over the hull form in 
the region of a sonar dome. While there was some apprehension 
about the criticality of changes to the geometry of the sonar 
dome, a telephone call to the Naval Sea Systems Command in 
Washington, D.C. {21], indicated that because of acoustic 
and hydrodynamic considerations the dome design should be 
maintained unchanged. 

The resulting procedure for modifying hull forms does 
provide good results for that portion of the ship below the 
design waterline. However, as outlined in chapter two and 
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Appendix C, there are situations where the method does not 
provide exact results, e.g., when the station of maximum 

beam and section of maximum area do not coincide. Another 
unresolved weakness of the modification scheme is that it 
still does not provide the degree of control over specific 
hull regions often desired, e.g., an attempt to preserve 

the configuration of a sonar dome will result in preservation 
of the centerline profile only, the three dimensional geometry 
of the dome will be uncontrollably altered. 

In summary fee has been concluded that the modification 
technique holds a great deal of promise for use with automated 
methods. In particular, the procedure as it currently exists, 
will provide excellent results when dealing with ships for 
which there is no rigid requirement to keep a specific region 
fixed. Not only are the desired coefficients and characteris- 
tics obtained, the resulting hull forms appear to be acceptably 
fair. 

As with virtually all work of this type there is still 
need for additional development. Specifically, it is felt 


that those aspects worthy of attention are: 


1. Investigate a means of controlling the resulting 
hull form above the design waterline. At present, 


excessive flair or tumblehome frequently occurs. 
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Investigate a means of rigidly controlling the 
geometry of a specific region of the hull. fMThis 
would provide a solution to the problem of 
keeping the sonar dome unaltered. 

Develop a computer program to carry out the 
extensive mathematical and graphical calculations 


required by the method. 
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6.2 Mathematical Representation of Lines and Fairing 


6.2.1 Lines Representation 


In chapter three it was demonstrated, by variational 
calculus and by simple beam theory, that a third degree or 
cubic polynomial could be used to approximate the shape taken 
by the draftsman's spline. However, it was also pointed out 
that the simple cubic polynomial became indeterminate if the 
curves contained infinite slopes. For this reason, and also 
because they are capable of representing multivalued functions, 
the parametric cubic equations of reference [15] were incor- 
porated into this thesis. The results obtained using this 
method have proven to be excellent. Not only does the tech- 
nique lend itself readily to being programmed, the parametric 
form of the curve allows the user to define either variable 
as being the independent variable for the purposes of inter- 
polation. The benefit of this capability will become 
apparent in the discussion of cross fairing. 

The only disadvantage, as seen by this author, to 
using the parametric equations is that they require the user 
to calculate the, up to three real roots, of the polynomials 
each time an interpolated value is sought. However, this 
is not seen as being restrictive since a closed form solution 


exists for calculating these roots and is in fact utilized in 
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subroutine CALCT. Therefore, because of its great flexibility, 
the parametric, or rotating spline technique of chapter three, 
is highly recommended for use with a lines fairing scheme 


involving the manipulation of specific waterlines and sections. 


6.2.2 Fairing 


The least-squares fairing criteria, as presented in 
chapter four, has shown to provide an effective means of 
altering the position of data points in order to obtain the 
desired "fairing" effect. That is, if provided with an 
adequate tolerance interval, the cubic spline passed through 
the resulting points will be void of extraneous oscillations 
and generally pleasing to the eye. When addressing the lines 
of a ship in the preliminary design phase, the fact that the 
lines satisfy a visual inspection is likely to be sufficient. 
For this reason, and also the excellent results obtained by 
this method in reference [16], this author has concluded 
that this scheme would be a candidate for a complete lines 


fairing program for destroyer-type ships. 


6.2.3 Recommendations 


It is obvious that, given the capability of repre- 
senting lines mathematically and also a means to fair the 


points on a line, the next step would be to generate a method 
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which would fair, in the three-dimensional sense, and display 
an entire ship. This author has spent a great deal of time 
attempting to extend the methods of reference [16] in a more 
general form to accommodate the peculiarities which arrise 
in addressing displacement-type ships. The difficulty arrose 
from two sources. First, an attempt was made to treat the 
entire ship, i.e., the bow and stern were not truncated as 
was the case with other methods examined. Second, in trying 
to treat a large variety of ships, conveniently called dis- 
placement-type, the author was confronted with the problem 
of attempting to describe the myriad of lines of discontinuity 
which one may encounter. These lines are most frequently 
termed control lines and may consist of the ship's profile, in 
an obvious sense, to the locust of points, longitudinally, 
where rise of floor and bilge radius meet, in a more subtle 
sense. Figure 6.1, for a typical bulbous bow destroyer 
illustrates a few of the possibilities. 

If we were to ignore the fairing algorithm itself 
for a moment, it can be seen that if a waterline A-A is taken 
in figure 6.1A there must be some means of communicating the 
effect of control line #6 on the waterline; where the explaina- 
tion of the control lines is contained in table 6.1. For this 
case, the effect is to create a straight line region in A-A 


as projected in figure 6.1B. A tentative solution to this 
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TABLE 6.1 


Explaination of Control Lines 


Bow profile 

Locus of stern radius centers 

Sonar dome profile 

Main deck centerline profile 

Deck edge profile 

Extent of deadrise 

Forward extent of parallel middlebody 
After extent of parallel middlebody 
Keelrise aft 

Outboard transom profile 

Transom centerline profile 

Deck edge waterline 

Forward extent of parallel middlebody 
After extent of parallel middlebody 
Extent of deadrise 

Outboard transom profile 

Deckedge transom profile 


Section view of transom 
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problem would be to ascribe to each control line, over a 
region where applicable, a code designating the effect of 
the control line on waterlines or sections at the point of 
intersection. This could be easily done by assigning another 
column to the CRV matrix description of the line, see section 
Sele3. 

Another complication which must be resolved is: 
when attempting to establish the offsets for, say an arbitrary 
waterline, how do you seek out where this waterline intersects 
which control lines. In the most general case, where control 
lines could occur at random through a hull form this problem 
could prove to be formidable at least. As seen by this 
author, the only solution to this problem is to have only 
certain control lines admissible for a particular class of 
ship. This would necessairly limit the possible intersection 
combinations. The control lines shown in figure 6.1 represent, 
what this author feels, are typical of a contemporary destroyer. 

The final aspect to be addressed is that of the cross 
fairing algorithm itself. Reference [16] showed that by 
utilizing a preassigned grid in the X-2 plane the offsets 
(y-coordinates) at these points could be repeatedly by faired 
and splined by both lines of section and waterlines. The 
new, or faired value of each point was taken to be the mean 
of that obtained by fairing the two lines. These mean values 


were then used as unfair data points on the lines once again 
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and the fairing process was repeated. This iterative pro- 
cedure was continued until the movement of the points on 
successive iterations was less than some predefined limit. 
The results of this cross fairing algorithm [16] proved to 
be quite good. Because of this, it is felt that this pro- 
cedure would also prove satisfactory for the more general 
method of lines fairing and representation presented in this 
thesis. 

As a final note, this author can envision where 
the two independent aspects of this thesis could be combined 
into one program of significant value. If both the fairing 
procedure and lines modification techniques were automated, 
it would provide the designer with the capability to sketch 
out a rough design on the back of an envelope, specifying 
its fundamental coefficients and dimensions, and then by 
passing this information through the fairing and modification 
routines a faired form could be obtained. The implications 
of this, as a savings of time and resources, are quite 
astounding. If the method were further extended to permit 
an interactive modification of the design, an individual 
could literally sit down and design a faired vessel ina 


matter of hours instead of days. 
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APPENDIX A 


Calculation of coefficient ¢ and centroid of the sliver of 
added area. Referring to figure 2.4. 
Recall: 6x = cx(l-x) 


i iL 
c f x(l-x)dy 
0 


Nl 
~ 
Or 
“ 
Qu 

i 


oo 


1 5 a 
e{fsfxdy- sf x"dy} = c[$-26x] 
0 0 
prey sen. 86 
o72oxX 6 (1-2x) 


gx = —S?_ x(1-x) (A.1) 
$ (1-2) 


solving for centroid, h 


a 


6oeh i 6x (x+8X) dy 
0 


ih l 1 
= f x6x dy + = ff 6x dy 
0 2 0 


substituting for 6x 


1 2 1 
5 2a3 8x 2 eses 
och = —S2 ff (x*-x") dy+—9¥ ft (x*- 2x" +x") ay 
o(1l-2x) 0 2(1=2x) 0 
= 2 2 2 3 
l - 2x o (1 = 2x) 
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For 6¢ << 6 


— 2 
Floss 2x 8 
1-2x 
where: 
ft i 
oe 3 3 aa 4 
k" = 39 : x dy r- = ip ; x dy 
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(A.3) 


APPENDIX B 


Sectionar Ares 
Curve. Au, nor 


fect SS. 





Referring to figure B.1 above, all longitudinal dimensions 
are measured with respect to the point of maximum sectional 


area. Assume: 


6x ex (1-x) 


8¢ 


4 a ee x (1l-x) 
 (1-2x) 


For the new hull form 


aa Sa 
gL = 5 (Lolo, + 665) + Co, + 66.) (B.1) 
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—-, _ 1 Dees m2. = 
z= Le {Le (¢.%, + 5h.) Ly (o,%, + 6¢h,) (B.2) 
t 


Also for the original hull form 
6, = E{L.¢, + L.¢,} eC) 
t L“e‘e ey Ais 
We now have to apply equations (B.1) and (B.2) to obtain 
values of So, and So. in terms of the known quantities z 
and dye These quantities representing the desired values LCB 
and C5 for the derived form. 


Solving equation (B.1) for 5o,° 


o,L = Lode Be L5¢, at L,. (9, + 6¢.) 


it : 
1 {Lo Sue Cees 6¢)} (B.4) 


56, 


Substituting equation (B.4) into equation (B.2) and solving for 


od _. 


ES 


a iam. LUO he ee 
2°L de = Lo (oe%~ + $$ hn.) L, ($x, + 6¢h.) 
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expanding the right hand side, R.H.S. 
L.H.S. = L26 x, + LESbh, - Led ¥, - L286 hn, 
rearranging terms 
L78¢ h_ - 1266 -h_ = L26_x. - L26x. - ZIL76! 
eo er e ‘ee e 
substituting in L.H.S. 


2 2 
L.o¢h,. = Loh 


i re ” = 
z Bete Lote L.(9,.+66,) } = R.H.S. 


2 o ; 2 
L $$ nh, LA lo, + Lhe. +) be nent 


Ir aS ec ek is LLh, 59, - 


2 
$$, (Lh +LLh,) : L, (Lh,o, ‘ L yhoo, i heUd,) a 


Se ee a eee 
5¢.. 2 {Lo _thy op %, zl + Loh, (Loe et 
Eon rh Lh 
2G as 6 ls SS) 


Equation (B.5) may be substituted back into equation (B.4) 


to obtain a value for oo. 
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However, if the simplified form for h is not used 
he= £(5o_,59.) 
Derivation of 8b, from equation (1) 


a sl ae - < 
= ue {Loy Lote 15, Lo} (B.6) 


So 


i: 
substituting into equation (B.2) and rearranging 


= pees 


ee ey | oe 
eee Oe Te%e e 


2 


7 Lh, 5%, a Lh {Lop -L,o,-L,5¢.-bb,} 


2 
= 6o {Lh +L Lh, }-LLh_ {Lo} -L¢-L.o,} 


4 1 2 
ove a ay! {Lio )X, La_X,t2'L ve Lh, (Le¢,t 
Lh+LLA 
ee sree r 
Lo Lot) } (BT) 


The form of this equation is merely the transposition of 


subscripts by r-e-r of equation (B.5) for So: 
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APPENDIX C 


Figure C.1-Section Modification 





Definitions: 
a = area of section 
b = beam of section 
dad = draft of section 


Subscripts: 
parent hull form 


PaS 
d = derived hull form 
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x = a maximum value, e.g., Dax is the maximum value 
of the beam in the derived hull form 
Superscript: 
bar, (—): refers to the sections in the parent and 
derived form which interact to create the 
derived section of maximum area. 
Pear, (*): Pe rers to the sections which are used to 


create the section of maximum beam in the 


derived form. 


Further define the following ratios: 


A= =, B=2, p= > 
x x x 
; 7 
_ Ag mead eed 
Oe A B=5- A a 
Pp Pp p 
also 
Bi) c: 
Bp 


It is by selecting a section in the parent, whose 
value of R is equal to that of the derived section being 
sought, that the new sections are created. -It should also 
be evident by referring to figure C.1, that the area of 


the derived section will be as follows: 


ie 





It remains to be shown that in some cases the resulting area 
ratio is not always what is desired. 


The following expression will also be useful. 


For any derived section, the area ratio obtained is: 


a b d 
oe: pik a Sh { — ao (Can) 
ax P px px 2-7 ax - dx 
anf Bo aa 
px px 
= _Baieipd: 
= ao 
p B A 
However, the area ratio desired is: 
a 
a = bs ds a 
RBs = po eA (ea2) 
p dd b j Day day ax 
Dox dox 


Therefore, the ratio of Aa desired to Aa obtained is: 


a 
aCe _ ae 
a a Bb. D 
fe OE aor = ee aS 
A a fox B_ D. 
ES 4 Sox px BB 
= B 
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Define: Ss = 


wiih 


R BD 
ppd 


It can be readily seen that, if the sectional area 
curve and the design waterline have their maximum values at 
the same longitudinal position, the value of S will be l, 
i.e., S = 1. Hence, the area curve obtained will be equal 
to that desired. 

If this is not the case, the designer has one option 
which will permit him to create a ship with the desired values 
of Cyr LCB, C., and LCF. 

*The designer must select a common longitudinal 
position about which to alter both curves. 
This will permit him to freeze either the 
point of maximum section of the point of 
maximum beam. Not both. 

*The only other alternative is to carry out 

the original procedure and accept slightly 
different values of C5 and Cy LCB and LCF 
will be as desired. The factor by which Cy 


will differ is: 


* 
me 2 (c.5) 


ors 





APPENDIX D, [15] 


Figure D.t - Rotating 
Spline Routine 





For n intervals bounded by n+l points the curve 


between points P. and Peal may be computed as follows: 


iewecemputce initially: 


arc tan [(yo-¥,)/(%5-x)) 1, if x) # Xo 


N/ 2; alee x, = Xo 


be P) = Jy 
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2. Compute n times for i = 1(1l)n 


z Bean weer 
Qe 8, HE(xy yp — 8) + Yay 7 YQ) 


iter ea) OG Mian) iV ig) ay) 





only if i> l 


Wet i =) 1 


a, on ee i et 


ce, 2) 27, 


ae 3k; (pj-9,) 
az k +2 
p,-r, (1+1/k;), ae k. x 0 
Be 
39;-Pi, 
<a! if k = 0 
3. Compute once: 
a. Gna. = 0 
b. d 


n+1 ~ Phtl a Gn+l 
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4. Compute n times for i = 1(-1)n 


ram ej eerste y 


a 
» 
I 
ctr 
y 
re 
a 
- 
t 
wQ 
- 


d. b.: 


Fi tan(d; 41-94) 


1 


The above procedure applies to the case where the ends 


2 


are pinned, i.e., d*y/dax = 0 at x = x, and x = x eee 


n+l 


however, it is desired to have the beginning slope equal to ty 


the following changes must be made: 


l 
ct 


egn. 1.b. P) = 


If it is desired to specify the end slope as ta41? 


See eyes tenes Pea 


rie 





To interpolate any point on the curve the following 


parametric equations are used: 


* 
I 


= X,t(%4 447%) t= (¥ 4g) -¥,) €(1-t) La, (1-t) -b, £1 
Y = ¥ytyy,17¥g) t= (%y477%;) t (L-t) (a, (1-t) -b, t] 


for 0=t<=t 


Section 3.2.4 describes the actual interpolation procedure. 
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APPENDIX E 


The following procedure was taken from reference [20] 
and is that used to obtain the roots of the parametric 
equations for x and y shown in Appendix D. 

Given the general form of the cubic polynomial: 


oe eee + bx + c= 0 (E.1) 


this may be reduced to the following by dividing by 


x =X) 7 a/3: 


x) = Ax, + B (E.2) 
where 

A = 3(a/3)7 -b 

B = -2(a/3)? + b(a/3) - ¢ (E.3) 
pefining 

p = A/3 and q = B/2 (E.4) 
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The roots of equation (E.2) are as follows: 


Case I: q? = p? > 0, there is one real root 


1/3 1/3 
= {q + ¥q*-p?} puna = Warep?) (E.5) 


ma 
There are also two complex conjugate roots. 


Case ITI: q? = p? = 0, there are three real roots of which 


two are repeated, i.e., only two roots are unique. 


am v3. ie 3. a 

xy = aa) a a) oad (q) v %4 a Xo (E.G) 
Case III: q? ~ p? < 0, there are three real and distinct 
roots. 

x, = 2¥p cos (U/3) 

Be 2¥p cos (U/3 + 2m/3) (E.7) 

ae 2¥p cos (U/3 + 41/3) 
where 

cos U = g/prvp 

olusn (E.8) 


is 





NOTE: These are roots of equation (E.2). To obtain the 
roots of equation (E.1) -p/3 must be added to the above 


solutions. 
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APPENDIX F 


Figures F.l and F.2 are conceptual flowcharts of 
the fairing and the splinning and interpolation procedures 
respectively. In the program listing that follows there is 
a short MAIN segment that requests the input data for a 
specific line and generates twenty-one (including end points) 
interpolated data points. While this program segment might 
prove of some value, it was designed primarily to test the 


various subroutines. 


NOTE: The program as listed requires the use of the LEQTIF 
Subroutine from the IMSL library. This subroutine is used 
in STRIP1 and STRIP2 to solve a 4x4 and a 3x3 system of 
simultaneous linear equations. For these SMAT is the 
coefficient matrix and T is the resultant column vector. 

If this library is not available any equivalent procedure 


could be substituted. 


NOTE: Due to time constraints at the time of publication, 
the program, as listed, will not accommodate curves with 
point types three or five. It will, however, handle curves 
without straight line segments and infinite slopes at end 


points. 
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CRV Matrix 


"DREFAR" 
Loads points on continuous 
segment into X, Y¥ and YORIG 


"FARCRV" NOTE; “FARCRV" 
Rotates axis to place X, calls "SPLINE" to 
Y and YORIG*XPRIM, YPRIM calculate the 
and YOPRIM Slope at point 3. 


"PARLIN" 
Loads or 5 points 
into X1, Y¥1 and YO 


? 
Are they 
irst points 


"estprs'" NO 
Calls STRIP1, 2 or 3 
to fair first 3 points 


Zz 
Are they 
ast points 


"TRANS1 " 
Transforms abscissa 
and calls STRIP1, 2 or 3 
to fair last 3 points 


Points are inte 
erior points 
and STRIP1 is called 


Figure F.1 — Flowchart of Fairing 
Subroutines 





124 





CRV Matrix 


"PRESPL" 
Loads points on a continuous 
segment into X, ¥ and YORIG 


Places A, B and D into the 
CRV matrix 


"SPLINE" 
Executes the rotating spline 
algorithm, calculating A, B and D 


"INTERP" 
Determines the interval 
in which a specified 
abcissa lies 


CALC Y 
Calculates the desired 
ordinate and slope at the 
specified abscissa based 
on the value(s) of 
the parametric variable T 


TOsneT: 
Calculates the real 
root(s) of the 
parametric variable T 


Figure F.2 —— Plowchart of Splinning 
and Interpolation Subroutines 
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APPENDIX G 


The following is an example of a typical bow section 
for a ship with a bulbous bow mounted sonar dome. Although 
the section is for a hypothetical ship it illustrates the 
capability of the computer program to fair, spline and inter- 
polate a curve having an infinite slope at an end point. 

In the CRV matrix shown, the items written in block 
numbers are input values while those in italics are values 
calculated by the program. The elements left blank were not 
used in this example. The values of tolerence, accuracy and 


limit are: 


TOL = 1.00 
ACC = 0.05 
LEiET = 10 
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+ Original Position 
e Faired Position 
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Figure G1-- Bow Section (Bulbous) 
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ao FF. WwW 


12.0 
17.0 
23.0 
28.0 
30.0 


CRV_ MATRIX 
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5 
G-3o3 
C2530 
0.394 

Ones 
-0.242 
= Oiae 106 
-0.041 
0.006 


6 
-0.480 
=e 
= eee 

0.160 
0.241 
0.149 
0.009 
-0.003 


eared. 
nine 
Ge 


- 366 
416 





The resulting interpolated values for increments of AX = (30-0)/ 


20 are as follows: 


ate DY/DX 
0.0 0.0000 Eee RA 
1.5 5.4206 1.3048 
3.0 6.8443 Orca oo 
4.5 7.4269 Oar ao 
6.0 7.3018 “052507 
7.5 6.4713 =O5 5019 
9.0 4.9988 =( 29816 
LO. 5 Saf 259 “007/295 
12.0 2.8047 -0.5114 
12-0 221714 =Oanoo 7d 
15. 0 1.7846 -0.1815 
16.5 Poot 9d -0.0415 
18.0 1.6542 0.0867 
19.5 bog72i 0.2010 
21.0 o3 egg 0.2956 
22.59 2.7465 0.3656 
24.0 3.335176 0.4116 
2on0 3.9724 0.4396 
27.0 4.6403 0.4476 
28.5 a. 9054 0.4368 
30.0 5.9543 0.4308 


Figure G.1 shows the curve generated as a result of fairing 


the given data points. 


159 


TOONS? 
Haberlandt 


Computer aided 


geometrical variation 
and fairing of ship 
hull forms. 





omputer aided geometrical variation and 


iit 





